Bridging the gap by shaking superfluid matter 
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In cold compact stars, Cooper pairing between fermions in dense matter leads to the formation 
of a gap in their excitation spectrum and typically exponentially suppresses transport properties. 
However, we show here that weak Urea reactions become strongly enhanced and approach their 
ungapped level when the star undergoes density oscillations of sufficiently large amplitude. We 
study both the neutrino emissivity and the bulk viscosity due to direct Urea processes in hadronic, 
hyperonic and quark matter and discuss different superfluid and superconducting pairing patterns. 



The dense matter in a compact star has certain trans- 
port properties, such as bulk viscosity and neutrino emis- 
sivity, that are dominated by beta (weak interaction) 
equilibration processes. However, this equilibration rate 
is expected to be exponentially suppressed as the star 
cools below the critical temperature for superconductiv- 
ity/superfluidity, and Cooper pairing produces a gap in 
the relevant fermion excitation spectrum jl] This 
suppression is relevant both in hadronic matter, where 
the critical temperature for proton superconductivity and 
neutron superfluidity is in the range Tc = 10^ — 10^" K 
[3] and in exotic phases such as quark matter, where crit- 
ical temperatures as high asTc^S — 5x lO^^Kif flavor 
antisymmetric pairing channels are favored , and as low 
as Tc ~ 10"^ K when flavor singlet pairing is favored [4^. 

In this letter we show that the exponential suppres- 
sion operates only at small amplitude and may be 
overcome in realistic situations. It is already known 
that beta-equilibration rates in normal matter are en- 
hanced by suprathermal processes at high amplitude 
when equilibration is driven by higher order terms in 
the chemical potential imbalance [5,-i7j. Here we report 
on a similar but more dramatic phenomenon in super- 
fluid/superconducting phases of matter, arising from a 
threshold-like behavior with a rapid increase in available 
phase space when the typical energy in the equilibration 
processes approaches the gap. 

The mechanism we discuss is generic and can be ex- 
pected to operate in all situations where large pertur- 
bations drive the system out of beta-equilibrium. One 
possible application are unstable oscillations of rotating 
compact stars, like f- or r-modes [8j. These modes are 
unstable due to gravitational wave emission and their 
amplitude grows exponentially until they are saturated 
by non-linear coupling to damped daughter modes [9 or 
damped by the supra-thermal bulk viscosity [10]. Note, 
that even at amplitudes well below the level where bulk 
viscosity alone would saturate such modes, suprather- 
mal effects can noticeably affect the thermal evolution 
of the star, via viscous heating and suprathermally en- 
hanced neutrino emissivity [6^ . Another class of scenarios 
stems from large amplitude oscillations caused by singu- 
lar events, like star quakes [Til or tidal forces preceding 
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Figure 1: Schematic illustration of the opening of phase space 
for Urea interactions in gapped nuclear matter when the de- 
viation from beta equilibrium /ia = Mn — A'p~Me becomes suf- 
ficiently large relative to the pairing gaps A„,Ap. 

neutron star mergers |12| . In this case the suprather- 
mal bulk viscosity should control the amplitude of the 
oscillation and correspondingly the emitted radiation. 

Here we will study the effect of large amplitude os- 
cillations in the case of various gapped phases of dense 
matter. When, as in nuclear matter. Cooper pairs con- 
tain two particles of the same flavor, density oscillations 
cannot lead to pair breaking. (Thermal pair-breaking 
has been used [T^ [TJ to explain the rapid cooling ob- 
served in Cas A [15].) Nevertheless, oscillations displace 
the gapped Fermi seas, and at sufficiently large ampli- 
tude can bridge the gap, opening phase space for beta 
equilibration processes as illustrated in flg. [ij As we will 
show, this suprathermal effect, which begins at ampli- 
tudes far below the level where non-linear hydrodynamic 
effects, of higher order in the fluid velocity, would arise 
[51|7], causes a huge increase in the corresponding rates, 
effectively bringing them up to their ungapped level. 

Beta equilibration processes have a finite time scale, 
so a harmonic local density fluctuation with amplitude 
An around the equilibrium value h induces an oscillation 
in the displacement /iA of the chemical potentials of the 
degenerate particles from their equilibrium. This oscil- 
lation leads both to heat generation (via dissipation due 
to bulk viscosity) and heat loss (from enhanced neutrino 
emission). Which of the two effects dominates depends 
among other factors on the amplitude of the oscillation 
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[B]. In [7] we found that over the entire range of physically 
reasonable density amplitudes the chemical potential os- 
cillation is given by the linear, harmonic relation 
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where C is a susceptibility characterizing the particular 
form of dense matter. 

We are interested in the direct "Urea" processes d — > 
u + I + Di and u + I ^ d -\- i>i mediated by charged cur- 
rent interactions, which occur in hadronic matter at high 
densities where proton fraction x > 10% |16^ To study 
all relevant cases within a unified framework we use a 
notation where d represents either an elementary nega- 
tively charged quark (down or strange) or a hadron that 
contains that quark. Similarly u stands either for an ele- 
mentary up quark or a hadron containing it; I represents 
a charged lepton (an electron, or in dense neutron mat- 
ter a muon) and h'l is the corresponding neutrino. We 
specialize to Fermi liquid theory, where the general dis- 
persion relation of the matter fields is 
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where i labels the matter species and represents a 
gap in the particle spectrum arising from superfluidity 
or (color)-superconductivity. Non-Fermi liquid effects 
should not play a role in a gapped system [T7]. We 
neglect dependence of the gap on energy and momen- 
tum, but we include temperature and density depen- 
dence Ai — Ai{T,fj,). The leptons can to very good 
approximation be described by a free dispersion relation 
Ef —pf + mf, =Pv Beta equilibrium with respect 
to the Urea processes enforces //a = ^J-d — fJ-u — fJ-i — 0- 
Since compact stars contain ultradegenerate matter with 
T ^ /i we can restrict the analysis to leading order in 
T/n where we find the general results for the net rate 
and the emissivity 
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or masses opens phase space and yields a non-zero re- 
sult. Eqs. ([3]) and ^ also contain a threshold function 
Q ~ 0{pFu+PFi—PFd)i a characteristic temperature and 
amplitude dependence, and a modification function R 
given for singlet gaps by the dimensionless integrals 
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where TVp : 
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functions reflect the phase space available to Urea reac- 
tions. They are normalized so that i?(0, 0, 0) = 1 and are 
symmetric in the gap parameters R{x,y,z) =R{x,z,y). 
In equilibrium they are pure reduction factors i?(0, ■ • • ) < 
1 that can become extremely small [T]. However, out of 
equilibrium ^a 7^ these modiflcation functions exhibit 
a suprathermal enhancement for ^a/T> 1 |[5jiZj and can 
then greatly exceed one. They depend on the equation of 
state only via dimensionless ratios and are the same for 
all direct Urea processes [Ij . In the ungapped case they 
have an analytic polynomial form [B] . 

The oscillation period of a compact star is much 
smaller than its evolution time scale, so the rel- 
evant quantity is the averaged emissivity tdu = 
5^ Jq^ (if ^duifJ-Ai^)) which takes the same form eq. Q 
with the averaged modification function Rg depending on 
/t A • The bulk viscosity of dense matter is [7] 



C = -^^ f^dipcosi^p) /Vr(^HMA(<^')) (7) 

Inserting the rate eq. ([3| and the low-amplitude form of 
the chemical potential oscillation eq. ([ij gives 



. _ 17G2 r4 ^f,^ Ad A„^ 



In these expressions G is the effective coupling of the 
fields to the weak current, D \s a. function depending on 
the density and the equation of state 
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which is identical for both quantities, and the parenthe- 
sis vanishes for a free massless dispersion relation. Aside 
from the neglected temperature corrections any modi- 
fication of the dispersion relations due to interactions 
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In hyperon and quark matter the viscosity is dominated 
by neutral current rather than Urea processes, but for 
those we also expect a high-amplitude enhancement |7] . 
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Figure 2: The modification functions Ri eqs. as a function of oscillation amplitude p.A, for different pairing patterns 

with maximum gap A: symmetric case where only one particle is gapped, e.g. Ad~A, Au=0 (solid lines); intermediate cases 
Atj = 5A„ = A (long dashed lines); Ad = 2A^ = A (short dashed lines); symmetric case A£j = Aii = A (dotted lines). Gap ranges 
from A/r = to A/T=100. Left panel: Ri for the averaged neutrino emissivity. Right panel: Rc^ eq. (|9| for the bulk viscosity. 
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Table I: Parameters for direct Urea processes in dense matter, with Fermi constant Gf, Cabbibo angle 6c, vector and axial 
couplings fv and qa [I], proton fraction x, symmetry energy S{n) |16| . in-medium hadron masses m*(n) and quark interaction 
parameter c [71[TS]. The quark expression is to leading order in Os and all results to leading order in T/n, mi/fi and pa/Mb- 



To see the gap-bridging enhancement of beta- 
equiHbration rates for oscillations of suprathermal am- 
plitude, without specifying a particular form of matter, 
we show in fig. [2] the modification functions Rg and i?^. 
We vary the maximum gap A and the ratio of the two 
gaps. For strongly gapped particles A^/T^l, the larger 
gap tends to dominate, e.g. (A^, A„) — (2A,0) is more 
suppressed than (A, A). Once the amplitude fi^ be- 
comes comparable to A, the Urea reactions can bridge 
the gap(s) so that the modification functions rise steeply, 
and quickly reach their ungapped levels. Interestingly, for 
matter where the two gaps are the same (dotted lines) 
this rise starts already for fi/^ <C and eventually merges 
into the solid curve for a single gap Ad+A.^- In contrast, 
in the asymmetric case where one gap is much larger than 
the other, the effects of the smaller gap are only visible at 
large amplitudes and become negligible when it is more 
than an order of magnitude smaller. 



Fig.|2]shows that suprathermal enhancement in gapped 
matter is considerably bigger than in ungapped matter 
(the A/T — line) [7 . This enhancement could be re- 
alistically achieved in compact star oscillations: fi/^ is 
related to the density amplitude An/n by the suscepti- 
bility C eq. (IT|) . For hadronic matter with an APR equa- 
tion of state [19] , C increases from C (no /4) w 20 MeV 
to C (Stiq) « 150 MeV P]; for quark matter it decreases 
from C(no) « 30 MeV to C(5no) ~ 20 MeV. Thus for 
amplitudes An/h ^ 0.01 the amplitude /tA can indeed 
become large enough to bridge gaps A ~ 1 MeV • 

To apply these general results to explicit processes re- 
quires dispersion relations and effective couplings for the 
relevant particles. Here we consider direct Urea processes 
in hadronic [2l[T6l|20], hyperonic [21121] and quark matter 
[221 [53] . Parameter values are given in tab. |l] Inserted in 
the general results eqs. ^ and ([8| they reproduce pre- 
vious expressions in the subthermal case fi^ <C T. In the 
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Figure 3: Enhancement via suprathermal oscillations of the 
neutrino emissivity and the bulk viscosity of hadronic matter 
with a ^5*0 proton gap Apo ~ IMeV at n « 1.3 no and a 
neutron gap A„o~0.1MeV at n~4no. 



gap at lower density is bridged, boosting the results by up 
to 10^°! As in previous results for ungapped matter [7], 
these effects are expected to be even more pronounced for 
the realistic case of modified Urea reactions [TJ [B] , that 
are actually present at moderate densities. 

Fig. [3] shows that, if the star is cold enough that Urea 
processes are blocked by pairing gaps throughout the 
star, density oscillations could provide the dominant con- 
tribution to those processes, by bridging the pairing gaps, 
starting in the regions where the rate-controlling gap is 
the smallest. The smaller this minimum gap, the lower 
the amplitude required to bridge it. In hyperonic pi I21| 
or quark matter |25l EB] there are processes which are 
only suppressed by A < 0.01 MeV, that could be bridged 
by oscillations with amplitude as small as An/n<lG^^. 
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support under contracts #DE-FG02-91ER40628, #DE- 
FG02-05ER41375, #DE-FG02-00ER41132 and by the 
Topical Collaboration "Neutrinos and nucleosynthesis in 
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ideal gas approximation our results simplify to 



TSRc (10) 



where Tg is the temperature in units of 10^ K and the 
values of e and ^ are given in tab. |l] 

As an explicit example that takes into account the com- 
plications in gapped interacting matter we consider APR 
hadronic matter |19| with the in-medium hadron masses 
given in For the density dependence of the ^5*0 pro- 
ton gap we use a Gaussian fit to data from [3], with 
maximum Apo — 1 MeV centered at n « 1.3 rig- For il- 
lustrative purposes we assume the neutrons also have a 
^Sq gap (in reality it is expected to be "^^2) with the 
density dependence proposed to explain the Gas A cool- 
ing data in [TJ, i.e. with maximum A„o = 0.12 MeV at 
n sa 3.7 riQ. The gaps have a BGS temperature depen- 
dence, but this has only a modest effect in small regions 
where r~ Ap,A„. Our results for neutrino emissivity 
and bulk viscosity at T = 10^ K are shown in fig. |3] At 
infinitesimal amplitudes (thick An/n— lines) there is 
enormous suppression by the proton gap, and moderate 
suppression by the smaller neutron gap. As the density 
oscillation amplitude increases, it first, at An/n > 10"'^, 
overcomes the suppression of emissivity and viscosity due 
to the neutron gap, increasing them at densities poten- 
tially relevant for direct Urea processes by many orders of 
magnitude. Then, at An/h>10~^ , even the large proton 
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